Resonant tunnelling between Luttinger liquids: solvable case 
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We discuss the conductance of a Luttinger liquid interrupted by a quantum dot containing a 
single resonant level. Using bosonisation and re-fermionisation methods, we find a mapping to 
a Kondo-type problem which possesses a non-trivial Toulouse-type solvable point. At this point, 
we obtain an analytic expression for the non-linear current-voltage characteristics and analyse the 
differential conductance and the width of the resonance peak as functions of bias and gate voltages, 
temperature, and barrier asymmetry. We also determine the exact scaling function for the linear 
conductance. 
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The field of one-dimensional interacting metallic sys- 
tems recently experienced another revival as single-wall 
carbon nanotubes (SWNTs) have been found to display 
transport properties consistent with the Luttinger liquid 
(LL) theory M. While the electrical transport through 
clean SWNTs has been investigated in different indepen- 
dent experiments, the transport properties of SWNTs 
with impurities (or in more complicated set-ups) are still 
to be studied in detail. Progress was recently made in 
this direction: in the manufacture of quantum dots 
on the nanotube basis was reported. Surprisingly, the 
authors found that the transport is dominated by a co- 
herent transmission (or resonant tunnelling) in a wide 
parameter range. 

The presence of a resonant level is known to enhance 
conductance. Indeed, for non-interacting electrons the 
local level hybridises with the conduction band causing 
a Lorentzian shaped peak in the density of states (the 
conductance being related to the Breit-Wigner scatter- 
ing cross-section via the Landauer formula). Unless the 
system is exactly at resonance, the picture remains qual- 
itatively the same for the case of interacting electrons, 
even LLs ||. Therefore one expects the conductance to 
increase upon lowering the temperature. On the other 
hand, at low temperatures, the conductance is known 
to vanish (unless exactly at resonance) due to the effec- 
tive enhancement of backscattering processes specific for 
the LLs. Hence a non-monotonic behaviour of the linear 
conductance as a function of temperature. The limiting 
cases has been thoroughly studied in Refs. || [|, ||, |(| 
but the full description of the cross-over remains an open 
problem. 

Recently, Nazarov and Glazman (NG) calculated the 
cross-over conductance in the weak electron-electron in- 
teraction limit (when the LL parameter g is close to I) 
by using the Landauer type approach supplemented by 
renormalisation group Q. In this paper, we wish to dis- 
cuss the opposite limit of strong interactions, which is 
relevant for such systems as SWNTs. We shall concen- 
trate on the special value of the coupling g — 1/2 and 
present an explicit solution of the problem at this point. 
At this particular value of the LL parameter and when 



the resonant level energy is tuned to match the equilib- 
rium chemical potentials in the leads, the resonant tun- 
nelling process is marginally relevant and its amplitude 
increases logarithmically upon lowering the energy scale 
H . It turns out that at low temperatures in the linear 
regime (i.e. in the limit of small bias voltage) the se- 
quential tunnelling dominates the transport for g < 1/2, 
while above that value the resonant transmission wins 
over J3|, 0]. Hence, apart from being supplementary to 
NG results, the exact solution at g — 1/2 yields insights 
into the interplay between these two transport mecha- 
nisms. 

We model the system by a resonant level (which can 
also be regarded as a single state quantum dot, so we 
use both terms) coupled to interacting leads, which is 
described by the following Hamiltonian (we ignore the 
spin degree of freedom throughout the paper), 



H — Hk + Ht + He , 



(1) 



where Hk is the kinetic part, Hk — Aeftd + 
Yli=R l HobPili describing the the electronic degrees of 
freedom in the leads Ho[ipi], and the resonant level with 
energy A with the corresponding electron operators be- 
ing eft , el. The dot can be populated from either of the 
two leads (i = R, L) via electron tunnelling with ampli- 
tudes 7,, H t = S i 7i[d^V ; j(0) + h.c.]. Here He describes 
the electrostatic Coulomb interaction between the leads 
and the dot, H c = \ c eftel J^i i>\ (O)^i(O) ■ This inter- 
action is a new ingredient we have introduced, absent 
in j| and Q. It does not, however, affect the univer- 
sality as we shall show. The contacting electrodes are 
supposed to be one-dimensional half-infinite electron sys- 
tems. We model them by chiral fermions living in an 
infinite system: the negative half-axis then describes the 
particles moving towards the boundary, while the posi- 
tive half-axis carries electrons moving away from the end 
of the system. In the bosonic representation Hq [ipi] are 
diagonal even in presence of interactions (for a recent re- 
view see e.g. ||; we set the renormalised Fermi velocity 
v = Vp/g = 1, the bare velocity being vp): Ho[ipi] = 
(47r) _1 J elx [d x (pi(x)] 2 . Here the phase fields 4>i{x) de- 
scribe the slow varying spatial component of the electron 
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density (plasmons), ipl(x)ipi(x) = d x (f>i(x)/2 / jr y /g. The 
electron field operator at the boundary is given by[fL4|, 
^(0) = e^^^/^/v^iao, where do is the lattice con- 
stant of the underlying lattice model. Here g is the con- 
ventional LL parameter (coupling constant) || |). In the 
chiral formulation the bias voltage amounts to a differ- 
ence in the densities of the incoming particles in both 
channels far away from the constriction || . The current 
is then proportional to the difference betwee the densities 
of incoming and outgoing particles within each channel. 

To the best of our knowledge, Hamiltonian (|l]) can- 
not be solved exactly even in the g = 1 case as long as 
Ac remains finite. However, after a transformation of S 
and d operators to the spin representation of the form 
S x = (dt + d)/2, S y = - d)/2, S z = dU - 1/2, 

one immediately observes that the Ac term is analogous 
to the 5 z -spin density coupling in the Kondo problem. 
The latter is known to be explicitly solvable at a particu- 
lar value of the longitudinal coupling: the Toulouse limit 
(see e.g. Q). Let us perform a similar calculation. As a 
first step we introduce new symmetric and antisymmetric 
fields 4>± = (<f>L ± 4>r)/V2, which still fulfill the bosonic 
commutation relations. Then we apply the transforma- 
tion H' = WHU with U = exp(iS z (j) + /y/2g~) which 
changes the kinetic and the Coulomb coupling parts of 
the full Hamiltonian to (we drop a constant contribution) 

H' K + H' C = H K + (Ac/tt y% - v / 2/^)S z d x M0) , 
and the tunnelling part (terms containing 7^) to 



H' t = (2na ) 



-1/2 



1R& 



i(j>-/ s/2g\ 



S- , 



(2) 



where S± = S x ± iS y — d\d. At the point g = 1/2 one 
can re-fermionise the problem by defining new operators 



(3) 



which fulfill standard fermionic commutation rela- 
tions. With the help of the particle density operator 
V'iV'i = dx4>±/2ir we can immediately write down the 
refermionised Hamiltonian, 

H = Ho[iP±] + {\ c -2tt)2S z ^ + tP + + AS z 

+ S + (j L ->P- + iRipl) + (iLlpl + 7R^-)S- . (4) 

In the case of the symmetric coupling jl — 1r this 
Hamiltonian is similar to that of the two-channel Kondo 
problem and, at the Toulouse point Ac = 2-7T, can be 
solved exactly (out of equilibrium) using the method of 
Ref. |ll[] . The novel ingredient in the following analysis is 
the extension to the asymmetric case. To take advantage 
of the Toulouse point we set the Coulomb coupling am- 
plitude to 27r in what follows. This not only removes the 
four fermion interaction but decouples the '±' channels 
making the '+' channel free (i.e. decoupled from the dot 
variables) . 




FIG. 1: Linear differential conductance at the resonance 
A = as a function of temperature for different values of 
the asymmetry parameter a. 

As we already mentioned, due to the linear disper- 
sion relation, the current through the system is propor- 
tional to the difference between the densities of particles 
moving towards the dot and away from it in either of 
the channels. Due to the chiral geometry we then have 
I ~ 00 ) — '/'lV'lI 00 ) j which, being transformed 

to '±' channels, results in I ~ 00) — ij)_ijj_{oo). 

Since the '+' channel is free, it doesn't contribute to the 
above formula. As the ' — ' channel is also free when away 
from the dot, in order to calculate the current we only 
need to know the scattering matrix of '— ' fermions de- 
termined by Hamiltonian (Q). The chemical potential of 
the incoming particles is determined by the bias voltage. 
Hence, the current is given by (we measure voltage in 
energy units, i.e. set e = 1) 

I(V) = G J dujT{uj)[n F {LO-V)-n F {uj)] (5) 

where n F denotes the Fermi distribution function and 
D(u>) = 1 — T(u>) is the energy dependent penetration 
coefficient of the '— ' particles from x < to x > 0. The 
pre-factor Go = e 2 /h is fixed by the requirement that at 
zero transmission D(u>) = (or perfect transmission of 
the whole structure) one obtains the correct conductance. 

The easiest way to obtain the transmission coefficient 
is the equations of motion method. Since we have two 
types of operators: for the electrons of the '— ' channel 
and for the resonance level (we go back to the original 
<V , d operators), we need two equations of motion, 

id t ip-(x) = -id x ip-(x)+6(x){~f L d-"f R d il ), 

id t d = Ad + 7i V-(0) + 7^1(0). (6) 

Integrating the first one around x = 0we obtain 

#-(0 + )-V-(0-)] =7Ld-7fld f . (7) 
Acting with <9 2 + A 2 on both sides of this relation yields 
(d i 2 + A 2 )ftM0+)-V>-(0-)] (8) 
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■ lR )d t + zA( 7 2 - 7 2 (0) - 2 lRlL d t ^l (0) . 



Now we can insert into this relation the momentum de- 
composition of the field operator -0_ 



ik(t-x) 

2tt 



b k 



for x < 
for x > 



(9) 



Because the dispersion relation is linear, lo = vk = k, 
we can use cj as the momentum variable as well as the 
energy variable. Inserting Eq.@ into Eq.(^j) and using 
ip-{0) = [^-(0+) + V-(0")]/2 results in 



E(b u 



; ) = -il3+{ auJ + b u ) + iiial^ + blj , (10) 



where we introduced the following objects: E = A 2 — uj 2 , 
P± = [(1 - 2a)A ± w]/2, 7 = wy^l - a), and a = 
7l/(7l + 7b) (^ ne asymmetry parameter). From now on 
ui, A, the bias voltage V, and the temperature T are all 
measured in units of T = 7^ + j R . Considering in addi- 
tion to Eq.([l0|) its complex conjugate for — u> we establish 
a relation between the amplitudes of the incoming (a w ) 
and transmitted (b u ) particle fluxes. The transmission 
coefficient can then be read off as follows: 



T(co) 



(11) 



4 7 2 £ 2 



{E 2 + /3 2 )(£ 2 + f3 2 _) + 2 7 2 (£ 2 + 0_0+) + 7 4 ' 

This equation, accompanied by Eq.(^), provides all in- 
formations about the transport properties of the sys- 
tem and is the central result of this paper. The experi- 
mentally relevant quantity is the differential conductance 
G = dl/dV. At zero temperature, Eq.(||) considerably 
simplifies and, differentiating with respect to the bias 
voltage, one immediately finds that G/Go = T(V). In 
the case when the couplings between the dot and the 
leads are perfectly symmetric and one of the chemical 
potentials matches A, G reaches the maximal value of 
Go • This is a typical signature of the resonant tunnelling 
effect usually encountered in transport phenomena in 
double-barrier structures In fact, our system is a 

model for such a structure with one single state between 
the barriers. 

The interplay between the LL's enhancement of the 
backscattering at low temperatures (resulting in decreas- 
ing conductance) and the more standard Breit-Wigncr 
physics emerging in the resonant tunnelling can be seen 
in Fig.Q. As predicted in Refs.||, Q], in the symmetric 
case a = 0.5 and for A = 0, the conductance saturates 
at low temperatures to its maximal value. In the presence 
of an asymmetry G does not saturate any more and van- 
ishes as a power-law towards T = with the exponent 2. 
This value is equal to twice the density of states exponent 
v of the LL with an open boundary: v = 1/g—l || , which 
in our case is equal to 1. This fact indicates that in this 
regime the electrons are transferred through the system 
in a single stage process Q, so that the internal struc- 
ture of the dot does not matter any more. Contrary to 
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FIG. 2: The width of the resonant conductance peak as a 
function of temperature for different values of the asymmetry 
parameter a. 

Eq.(9) of NG, the high temperature (T ^> 1) evolution 
of the conductance follows the law G/Gq ~ \/T. The 
reason is that the problem maps onto a free-fermion one, 
for which the 1/T behaviour is inevitable. Note that, 
in the language of the original model, this corresponds 
to tunnelling of composite objects from one LL into the 
other. 

Another interesting issue is the shape of the resonance 
peak, especially its width, w(T), as a function of tem- 
perature, see FigJ|. At high temperatures it decreases 
linearly upon lowering T no matter how strong are the 
interactions. For T< 1, however, the correlation effects 
become visible and the width w(T) of the peak saturates 
at zero temperature unless the dot is symmetric. In the 
latter case w(T) shrinks to zero with the exponent 1 — g 
predicted in Ref.fJ: w(T) - T 05 . 

It is in fact not difficult to evaluate the integral in 
Eq.([5]) analytically. The general expression is compli- 
cated, so we shall only present here some particular cases. 
To begin with we observe that there is an intimate rela- 
tion between our model at A = 0, a = 1/2 (resonant 
and symmetric case) and the g — 1/2 solution for the 
conductance through a single barrier (G s ) given in ||. 
Indeed, evaluating ([ll]) for the case in question, we find 
(in linear response) 



(12) 



where tp is the ^-function. Comparing with Ref. p3[ (see 
also §), we observe that G A =o{T)/G = 1 - G S (T)/G 
if T in G s is measured in units of the backscattering 
strength. One can easily show that an analogous relation 
continue to hold for the out-of-equilibrium current. 

Furthermore, for the linear conductance when a = 1/2 
but A/0 (|A| < 1/2) we obtain 



G A (T)/G 
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FIG. 3: Comparison between the linear conductance (dashed 
lines) and the approximative scaling function (iu/A) 4 (solid 
lines) for different values of A: from above A = 0.5, 1, 2, 4. 
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(13) 



where \/2A ± = V 1 - 2A 2 ± \/l - 4A 2 . Same formula 
is valid in the resonant but asymmetric case (A = 0, 
a ^ 1/2) if we substitute A± -> 1/2 ± ^/a(l - a). This 
shows that the asymmetry parameter is equivalent to a 
small off-set of the resonance (if the two act indepen- 
dently). (Indeed, in terms of the Kondo analogy, the 
channel asymmetry is known to be relevant.) For a strong 
off-set of the resonance (|A| > 1/2), formula ( |i~3| ) is not 
valid and should be replaced by: 



G A (T)/G = 



7r7V4A 2 - 1 



-Am 



v(i + A ' 



J t) 



(14) 

where 9 = tan^^A 2 -1/(1 - 2A 2 )]/2. 

As pointed out in Ref.g, for A, T <C 1, the conduc- 
tance should become a universal (scaling) function of 
the ratio of the resonance width and the backscattering 



strength (resonance off-set). Indeed, taking the appro- 
priate limit in formula (fL3J), we obtain the exact scaling 



function at g = 1/2: G A (T)/G Q = G g=1/2 (X), where 



G 1/2 po = i--xVU + 



■A' 



(15) 



and the scaling variable is X = y / 7rA/2T 1 / 2 (that is 
w(T) ~ 2T 1 / 2 / v /tt at small T). We note that our scal- 
ing function G g=1 / 2 (X) is by far more complicated than 
G g ~i(X) = 1/(1 + X 2 ) found by NG in weak coupling. 
Furthermore we observe from our analytic expressions 
that beyond A,T <C 1 there is no exact scaling. How- 
ever, upon determining w(T) numerically and plotting 
[w(T)/A] 4 versus the dimensionless conductance G/Gq, 
we obtain an approximate numerical scaling as shown in 
Fig||. The same scaling function ( |l5| ) holds in the reso- 
nant (A = 0) but weakly asymmetric case, when a — 1/2 
is small. In that situation a — 1/2 substitutes A in the 
definition of X. 

To summarise, we presented an explicit solution for 
the transport through a resonant level coupled to two 
LL leads. It turns out that for g = 1/2 the Hamiltonian 
of the system can be mapped onto one similar to the 
two-channel Kondo Hamiltonian in the Toulouse limit, 
solvable exactly. We obtained the full I — V character- 
istics, which shows all the effects inherent to resonant 
tunnelling setups in LLs, including the scaling. Our so- 
lution confirms previous results obtained by means of the 
perturbation theory and goes beyond them. In future, it 
would be interestin g to study deviations from g = 1/2 
(in the spirit of Ref. |L3|) and to investigate the effects of 
electron spin (and flavour). 

The authors would like to thank H. Grabert and 
R. Egger for a valuable discussion. This work was 
supported by the Landesstiftung Baden-Wiirttemberg 
gGmbH (Germany), by the EPSRC of the UK under 
grants GR/N19359 and GR/R70309, and by the EC net- 
work DIENOW. 



[1] C. Dekker, Phys. Today 52, No. 5, 22 (1999), and refer- 
ences therein. 

[2] H. W. C. Postma, T. Teepen, Z. Yao, M. Grifoni, and 

C. Dekker, Science 293, 76 (2001). 
[3] C. L. Kane and M. P. A. Fisher, Phys. Rev. B 46, 15233 

(1992). 

[4] Yu. V. Nazarov and L. I. Glazman, cond-mat/0209090 
(2002). 

[5] A. Furusaki, Phys. Rev. B 57, 7141 (1998). 

[6] C. Nayak, M. Fisher, A. Ludwig, and H. Lin, Phys. Rev. 
B 59, 15694 (1999). 

[7] K. A. Matveev, Sov. Phys. JETP 72, 892 (1991). 

[8] A. O. Gogolin, A. A. Nersesyan, and A. M. Tsvelik, 
Bosonization and Strongly Correlated Systems (Cam- 
bridge University Press, 1998). 



[9] R. Egger and H. Grabert, Phys. Rev. B 58, 10761 (1998). 
[10] V. J. Emery and S. Kivelson, Phys. Rev. B 46, 10812 
(1992). 

[11] A. Schiller and S. Hershfield, Phys. Rev. B 58, 14978 
(1998). 

[12] S. Datta, Electronic Transport in Mesoscopic Systems 

(Cambridge University Press, 1995). 
[13] U. Weiss, R. Egger, and M. Sassetti, Phys. Rev. B 52, 

16707 (1995). 

[14] Strictly speaking ip(x = 0) = 0, so we assume that the 
tunnelling takes place at the second last site of the cor- 
responding lattice model, at x = ±ao. Also, we ignore 
the Klein factors as they can be absorbed into Eqs.(^) 
disappearing from the analysis. 



